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ABSTRACT: We consider a simple spin system without disorder which exhibits 
a glassy regime. We show that this model can be well approximated by a system with 
quenched disorder which is studied with the standard methods developped in spin 
glasses. We propose that the glass transition is a point where quenched disorder is self 
induced, a scenario for which the 'cavity' method might be particularly well suited. 



The problem of the glass state remains one major unsolved issues in condensed 
matter theory. Despite an enormous body of experimental and numerical data and 
quite detailed phenomenological theories [1,2,3], there is no fully satisfactory micro- 
scopic model for the glass state. The intense theoretical activity on spin-glasses and 
other disordered systems [4] stemed in part because they retain 'half of the complexity 
of glasses: given a disordered ('quenched') set of interactions, what is the thermody- 
namics of the 'spin' degrees of freedom, is there a low temperature spin glass phase, 
etc... The spin glass theory has indeed given birth to many seminal ideas which have 
been transfered to other glassy systems like proteins [5,6,7] , rubber [8], or even glass 
itself [3]. One subtle aspect of glasses is that there is no clear a priori distinction 
between 'slow' degrees of freedom responsible for random interactions and 'fast' de- 
grees of freedom equilibrating therein, although everything goes as if it was the case: 
quenched disorder is self-induced. A satisfactory glass theory requires a detailed math- 
ematical description of this scenario (in fact implicit in the mode coupling theory [9]) 
and the identification of these 'slow' degrees of freedom . 

In the following we shall show a possible way to get round this problem. We 
present a simple model, already studied in [10,11], which contains no quenched disor- 
der. At this stage, the only justification for studying this model comes from numerical 
simulations : the system very clearly exhibits features of a glass transition - jump in 
the specific heat [11], slow dynamics and aging [12], etc... We propose an unusual 
analytical approach to this model which is to find a 'fiduciary' disordered model which 
is 'as close as possible' to the pure model, but for which all the ideas and methods 
developed for spin glasses (replicas, cavity method, statistics of the metastable states) 
are readily available. Our approach is the complete opposite of the usual one, which 
is to replace a 'dirty' system by an equivalent, 'pure' one. We show that the high 
temperature (replica symmetric) phase of our 'fiduciary' system reproduces exactly an 
approximation due to Golay [10,11] for the original model, which, although unjustified, 
accounted reasonnably well for Bernasconi's numerical data at high enough tempera- 
tures [11]. The entropy given by this approximation however becomes negative at low 
temperatures, signalling, for our fiduciary model, the breaking of replica symmetry. 
We find that the system undergoes a first order transition towards a low temperature 
(glass) phase which is rather similar to the low temperature phase of Derrida's ran- 
dom energy model [13], although the entropy remains non zero - reflecting the fact 
that small scale motions are not completely frozen. This random energy structure is in 
good agreement with the numerical findings of Bernasconi, who found that the energy 
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landscape is 'golfcourse' like, with low energy states randomly distributed in phase 
space [11]. 

We shall first describe the specific model we considered and its fiduciary disordered 
version, and sketch the main steps of the calculations. We shall then turn to a more 
physical (and speculative) discussion on the relevance of the rather abstract model 
studied here for more realistic situations. 

The model in question is defined by the following Hamiltonian: 
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where Si = i _/y are Ising spins. The scaling of 7i with N has been chosen such that 
7i is extensive. The spin configurations which minimize 7i are binary sequences with 
small autocorrelations, which are useful in communication engineering problems [11]. 
It is difficult to find them because of frustration effects. 

As mentionned in the introduction, numerical studies show very clearly that the 
system enters a glassy phase at low temperatures, much as if quenched disorder was 
present. Furthermore, the non trivial features of H come from the fact that the sum 
over k extends to infinity when N — > oo. In other words, it is the couplings between 
very far away spins which matter - suggesting that the one dimensional nature of the 
problem might not be crucial. We thus propose to replace Eq. (1) by the following 
'fiduciary' Hamiltonian: 
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where JE are random connectivity matrices, independent for different /c's, with each 



element equal to Jq with probability 



N-k 



and zero otherwise. This choice insures that 



the average number of bonds in 7i and Ti^ is precisely the same: note that the choice 

(k) 

JE = Jo^i+fc ?' reproduces exactly Eq. (1). (Jq is set to 1 in the sequel). can be 

IJ'iJ 

considered as the mean field version of 7i where the geometry is lost. Interestingly, 
this mean field Hamiltonian allows one to use the replica formalism. After rather 
standard manipulations [4], we find that the free-energy at temperature 1//3 is 
given by — lim^^o where the average over the disorder of the n-th power of 
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the partition function is 
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Qa^b = <?• We find that g sad dle = Saddle = 0, leading to = -[5 1 J Q dx\og[l + 
(3(1 — x)]. Interestingly, this free energy coincides exactly with the one obtained by 
Golay [10] for the original model Eq. (1), under the (unjustified) assumption that 
Rk = Yli^i* SiSi+k are Gaussian independent variables. As shown numerically by 
Bernasconi, F^ s gives a rather good description of the 'high' temperature region. This 
solution however suffers from the usual entropy disease, which becomes negative below 
a certain temperature T* = 0.047564... and goes to — oo for T = 0. However, there is 
no sign of local instability, suggesting that the transition to a replica symmetry broken 
phase must be first-order (from the point of view of the order parameter function: as we 
shall see, the transition is second order from the ther mo dynamical point of view). The 
existence of a phase transition is ensured by the fact that Z 2 ~ at high enough T. 
The one step replica symmetry broken solution allows to introduce, as usual, a minimal 
and a maximal overlap qo,qi, as well as the position of the 'breakpoint' m, connected 
with the density of low-lying states (and in turn with the dynamical properties [14,15]). 
We find that q = and m(T) ~ ^ with T g = 0.047662 > T*, while y(T) = (3{l - q\) 
behaves as shown in Fig 1. Fig 2-a, 2-b show the free energy and the entropy in the 
low temperature phase T < T g . Note that the entropy is everywhere positive but 
rather small (~ 10 per spin), goes to zero linearly with T (as in real glasses), and 
matches that of at T = T g . T g thus appears as a freezing temperature at which 
gi discontinuously jumps from zero to a value rather close to 1 (Note the scale in Fig. 
1). The specific heat also jumps at T g . The picture of the glass phase is rather similar 
to that of the random energy model, for which [13,16] qo = 0, m(T) = if- but also 
q\ = 1 corresponding to the fact that the entropy of the frozen phase is strictly zero - 
at variance with our model for which a residual entropy remains. Our prediction for 
the ground state energy is Eqs = 0.02028455.... Bernasconi noted that the numerical 
ground state energy was much higher than this value as soon as iV > 50, which might 
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simply reflect the fact that if the energy landscape is that of the random energy model, 
it is extremely difficult to find the ground state. 

It could be that a more complicated replica symmetry breaking scheme is needed. 
We however think that the numerical difference with our results are likely to be ex- 
tremely small; it is furthermore irrelevant to the point addressed in this letter. We 
shall now discuss the above results from a more physical point of view. Although quite 
remote from reality, the model we considered illustrates the fact that a pure model can 
undergo a glass transition which can be described using the tools of disordered systems 
since its mean field formulation naturally introduces random variables. Breaking of 
replica symmetry indicates, as usual, the existence of many (metastable) states, and 
its physical meaning is best understood in the 'cavity' approach [4], which is essen- 
tially based on a certain (hierarchical) construction of the equilibrium states and the 
local field distribution. In disordered systems, this method is in fine equivalent to the 
replica calculation. In the pure model at hand, however, the cavity method in fact 
allows to recover all the above results without introducing a fiduciary random Hamil- 
tonian, but rather through adequate hypothesis on the statistics of the Si and the 
[17]. Within a one pure state picture ('replica symmetric'), the only viable assump- 
tion is that < Si >= 0, which immediately leads back to Golay's approximation. The 
existence of many 'states' a with weight P a however allows to go beyond this result, 
since it is possible to have ^ a P a < S^ > a = but J2 a < &i > a = 11 ^ 0- In other 
words, if the spins are frozen in a given state, the field acting on the extra ('cavity') 
spin will be much like a quenched random variable. A simplistic way to express this 
idea might be the following: one can rewrite Eq.(l) as 



with Jij = —^^kSi+fcSj+k, i.e. as a spin-glass SK Hamiltonian but for which 
the couplings are themselves determined by the spins. However, if the dynamics 
of the system is slow (which it is in the 'spin-glass' phase) then the couplings can 
be self-consistently thought of as quenched random variables. This suggests that 
the scenario described here is far more general [9] and that a genuine short range 
model of glass could be a pure four spin interation Hamiltonian of the form H4 = 
— J J2<ijkl> SiSj^kSh where < ijkl > denotes, for example, nearest neighbours tetra- 
hedra on a cubic lattice. H4 obviously posseses a 'crystalline' ground state S{ = 1. If, 
however, the system is quenched from high temperature, the {Si} are initially random 




(5) 
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and generate random effective couplings Jjj, which, if the temperature is sufficiently 
small, will very slowly evolve and lead to a 'self consistent' spin-glass. Numerical sim- 
ulations [18] show that this is indeed the case: glassy dynamics and aging very similar 
to that observed in experimental spin glasses is clearly observed. Of course, if the 
interaction is of finite range, this quenched effective disorder will progressively anneal 
out (possibly on astronomical time scales), allowing the system to find its crystalline 
state - as indeed in real glasses. Only if the interaction is of infinite range, like in 
the models considered above, or if some topological contraints forbid the annealing of 
disorder [19] will this glass transition acquire a precise ther mo dynamical meaning. It 
is clear, however, that the range of interactions need not be very large for this limit to 
be relevant to experimental time scales [20]: in this work, we have primarily focused 
on static calculations, leaving the investigation of the dynamics (along the lines of 
[21,22]) for future work. A dynamical approach is clearly needed in order to identify 
which internal degrees of freedom get quenched and to make a link with the mode cou- 
pling theory [9]. The success of the present approach might indicate that the precise 
decomposition of which degrees of freedom become quenched or annealed is maybe not 
so crucial. 
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Figure captions. 

Fig 1: Plot of the quantity y(T) = — q jS T ^ in the low temperature phase. Note 
the scale, which indicates that the deviations of q\ from 1 are very small. 

Fig 2: Plot of the free-energy (2-a) and entropy (2-b) in the glass phase. Note 
that S(T) is positive but rather small. We found that S(T) = AT (for small T) with 
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A ~ 5 x 10~ 5 
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